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MULTI-COMPONENT  MODEL  OF  THE  RESISTIVE  HOSE  INSTABILITY 

I.  Introduction 

The  resistive  hose  Instability  Is  a  growing  lateral  distortion  of  an 
energetic  self-pinched  beam  propagating  in  a  dense  resistive  plasna. 
iJhen  there  Is  no  significant  equilibrium  return  current,  the  instability 
Is  driven  hv  a  resistive  lag  of  the  magnetic  field  In  responding  to 
transverse  displacements  of  the  beam.  The  presence  of  return  current 
fiirther  destabilises  the  beam  by  magnetically  repelling  the  displaced 
beam  from  the  axis. 

Mthough  the  hose  Instahllltv  has  been  studied  extensively*  ,  most 
nrevlous  theoretical  worV  is  based  on  oversimplified  models  of  beam 
particle  dynamics.  ‘?lnce  the  pinch  force  Is  anharmonic  v^henever  the 
radial  profile  of  current  density  Is  rounded,  particles  have  a  spread  of 
betatron  frequencies  that  Introduces  phase  mixing  and  tends  to  damp  hose 
oscillations.  Tn  addition,  the  betatron  frequency  of  a  particle  depends 
on  the  extent  and  shape  of  Its  orbit,  leading  to  radially  localised  wave- 
particle  resonances.  Tlie  early  rigid  bean  hose  models^*  ^  and  I’.lnetlc 
treatments  assuming  helical  particle  orbits’’  omitted  phase  mixing 
entirely,  and  as  a  result,  these  models  overestimated  the  degree  of 
Instability.  The  "spread-mass”  model  of  Leo'’  was  a  significant 
improvement;  by  Introducing  a  realistic  amount  of  nbase  mixing  In  an 
artificial  way.  It  gave  a  reasonable  hose  growth  spectrum,  even  though 
the  model  neglected  the  correlation  betvmen  the  betatron  frequency  of 
particles  and  their  radial  location. 

l>ecentlv,  llhm  and  T^mpe^  developed  an  "energy-group"  model  that 
exhibits  radially  localized  resonances  as  well  as  nbase  mixing,  and  here 
we  present  an  Improved  and  refined  version  of  that  model.  As  In  the  Ilhm- 

Lampe  model,  each  beam  slice  Is  represented  by  a  superposition  of  rigid 
Manuscript  submitted  December  17, 1981. 
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components  with  different  outer  radii,  each  of  which  responds 
independently  to  the  resultant  transverse  force  actlni;  on  It.  Since  the 
hean  response  depends  sensitively  on  the  radial  density  profile  of 
components,  we  use  a  Vlasoy  analysis  here  to  ^lide  the  choice  of 
component  profiles.  With  an  appropriate  component  shape,  this  "multi- 
component"  beam  dynamics  model  accurately  reproduces  the  localized 
resonances  and  other  analytic  properties  of  a  full  Vlasov  treatment, 
while  still  beln^  tractable  by  straightforward  numerical  methods. 

The  outline  of  the  paper  Is  as  follows.  In  Sec.  TI,  r-re  formulate 
the  Vlasov  eigenvalue  problem  for  the  hose  Instability  and  determine  the 

A 

analytic  properties  of  the  perturbed  c»trrent  J{r).  The  details  are  given 
In  Appendices  A  and  B.  In  Sec.  Ill,  we  use  the  multi-component  model  to 
obtain  an  expression  for  J(r)  and  show  that  this  form  preserves  the 
important  analytic  properties  of  the  Vlasov  restilt.  The  method  of 
numerical  solution  for  the  dispersion  relation  and  eigenfunctions  Is  then 
outlined.  The  results  of  extensive  numerical  calculations  are  given  In 
Sec.  IV.  We  compare  hose  eigenfunctions  and  growth  rates  to  earlier 
models  and  discuss  the  sensltlvltv  of  the  growth  rates  to  the  equilibrium 
return  current  fraction  and  the  radial  profile  of  conductivity,  finally, 
we  mention  work  in  progress  that  extends  the  range  of  applicability  of 


II.  Hose  F.lgenvalue  Problem 

A.  Assumptions  and  Basic  Equations 

The  equilibrium  conslHereH  here  is  an  axlsymmetrlc  self-pinched  beam 

of  relativistic  electrons  moving  In  the  positive  z  direction,  with  no 

bulk  rotation  or  externally-imposed  fields.  We  make  the  "paraxial" 

assumption  that  the  transverse  velocity  of  beam  particles  Is  small 

compared  with  the  parallel  velocity  v^.  This  assumption  requires 

that  11  1  «  I.,  where  the  net  current  I„  is  the  sum  of  the  beam  current 
n  A  " 

> 

and  plasma  current  Ip,  and  I^  Is  the  Al^ven-Lawson  current,  given  In 
2  1/2  I  -1 

cgs  units  by  =  (y  -1)  Me  |q |  .  Here,  M  is  the  mass  of  beam 

particles,  q  Is  the  charge,  and  y  Is  the  typical  relativistic  factor.  It 

follows  that  the  beam  radius  a  is  small  compared  with  the  betatron 

1/2 

wavelength  of  beam  particles  X.  *  iTa(2I./|I  I)  .  The  background  gas  Is 

p  An 

taken  to  he  a  stationary  medium  characterized  by  a  conductivity 
oCr)  large  enough  to  provide  space  charge  neutrality,  which  requires 
that  c/(4xo)  «  a.  The  effects  of  collisions  between  bean  particles  and 
background  gas  particles  are  Ignored. 

In  the  hlgh-conductivlty  regime  considered  here,  a  non-rotating  beam 
has  only  an  axial  vector  potential  component  A^,  and  all  beam  fields  can 
he  derived  from  It: 


Rg  -  -  3A^/3r,  (la) 

B  -  r“*aA  /ae,  (ib) 

r  z 

E  ■  -c  ^aA  /at.  (ic) 

z  z 
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A  is  itself  determined  hy  Ampere's  equation. 


-7^  A  -  —  f  J.  J  ), 

i  z  c  h  p 


Here,  Jy,  is  the  7.  component  of  the  beam  current  density,  and  the  plasma 

current  density  .T_  in  (2)  is  Riven  hy  oE  •  Hor  a  paraxial  beam  in  this 
P  z 

regime,  transverse  currents  can  he  neglected  as  sources  in  the  field 
equations. 

The  eqnlllhrlum  distribution  function  f^,  can  be  defined  as  any 
function  of  the  constants  of  particle  motion,  the  particle  energy 


H  -  YMc^  i  (c^jj'R 


the  azlmithal  angular  momentum 


and  the  axial  canonical  momentum 


where  £  «  YMv  Is  the  particle  momentum  and  Af(  is  the  enulllbrlum  A^. 
Since  we  consider  a  non-rotating  beam  without  external  fields,  we  can 
then  take  to  be  independent  of  Pg.  To  model  a  paraxial,  nearly- 
monoenergetic  beam,  we  assume  that  all  particles  have  the  same  but 
there  is  a  spread  in  perpendicular  energy 
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5  H  -  y^Mc  , 


whare 


w  2  _  r  2  ^  A  a/2 
Y^Mc  =  (  Pj,c  +  M  c  J 


is  the  energy  associated  with  P  .  We  therefore  take  fy  to  have  the  form 


IqCH.P,)  -  «<P,-P,). 


It  then  follows  from  the  paraxial  approximation  that  v^  Is  nearly 

constant  for  all  beam  particles  and  approximately  equal 

to  0c  =  P  /(y,.M)»  For  this  distribution  function,  the  equilibrium  beam 
z  b 

current  density  is 


bW  f  ”  u(r:P„) 
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where  the  radial  velocity  v^  =  p^/ (YH)  is  obtained  as  a  function 
of  r,  H,,  and  Pg  from  (3)-(5),  and  the  ralnlraun  perpendicular 
energy  U(r;Pg)  in  (9)  is  calculated  from  the  requirement  that  v^  be 
real.  Since  all  beam  particles  in  this  model  have  nearly  tlie  same  axial 
velocity.  It  is  convenient  to  use  z  and  A  =  6ct-z  as  Independent 
variables  Instead  of  z  and  t.  Then  C  Is  a  constant  of  the  motion  and 
labels  a  co-movlng  slice  of  the  beam,  wlille  z  describes  the  time 
variation  In  the  particle  reference  frame. 

In  this  paper  we  are  primarily  Interested  In  developing  an 
eigenvalue  treatment  of  the  hose  Instability.  Since  this  approach  is 


possible  only  if  the  equilibrium  Is  Independent  of  C  and  z,  we  require 
that  the  beam  radius  a  and  the  equilibrium  beam  and  plasma  densities 
Jbo(r)  and  Independent  of  C  and  z.  These  quantities  are. 

In  fact,  slowly  varying  except  near  the  beam  head  *  ,  which  we  do  not 

consider  here.  In  addition,  we  must  ignore  the  C  and  z  dependences  of 
the  conductivity  o(r),  which  can  result  from  heating  and  Ionization  by 
the  beam.  ?luch  a  variations  can  have  a  significant  effect  when  the  beam¬ 
generated  conductivity  Is  a  large  fraction  of  a,  but  this  approximation 
is  appropriate  for  beams  propagating  in  highly  prelonlzed  channels.  We 
emphasize  also  that  the  multi-component  model  developed  In  Sec.  HI  Is 
not  limited  to  axially  uniform  beams.  A  multi-conponent  "simulation" 
code,  which  explicitly  retains  all  5  and  z  dependences,  has  been 

Q 

developed  and  will  be  reported  in  a  separate  paper. 

We  formulate  the  hose  eigenvalue  problem  by  decomposing 
perturbations  Into  Independent  Fourier  modes  with  azimuthal  mode  number 
m=l  and  axial  wave  number  k^.  In  terms  of  C  and  z,  the  perturbed  A^  and 

"^b  have  the  form 


(r,5,z) 

(r,?,z) 


exp(i9  -  1-^  C  -  i~  z). 


(10) 


where  SI  =  u>-k^0c  Is  the  frequency  seen  by  beam  particles,  and  a)/(Pc)  can 
be  thought  of  as  a  wave  number  for  spatial  oscillations  in  5.  From 
Ampere's  equation  (2),  the  perturbed  A^  satisfies  the  linearized  relation 


did  *  4Trlma(r)  ' 

T"  rA  + - -,7^-  A 

dr  r  dr  2 


(11) 


The  perturbed  beam  current  J  Is  obtained  from  a  standard  Vlasov  analysis, 
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first  Integrating  the  linearized  Vlasov  equation  along  unperturbed 
particle  orbits  and  then  Integrating  over  the  constants  of  the 
motion  and  Pg.  The  result®  Is 


2  2^  ^  A 

j(r)  -  /  dPg  /  dH^^  ^  [A(r)  +  M  i(r;H^,Pg)]  .  (12) 


where 


2 

I(r;Hj^,PQ)  -  I  I  /  dz'A(r;)  exp[l(0;  -  6)  -  l-g|(z'-z)].  (13) 


In  the  orbit  Integral  (13),  r_^  and  0^  specify  the  unperturbed 
location  at  z'<  z  of  a  particle  with  r^(z)  »  r,  0_j^(z)  “  ©, 

Vj.(z)  ■  +  |v^(z)I,  and  constants  of  motion  Pg  and 

B.  Particle  Orbit  Properties 

Particle  dynamics  In  the  transverse  plane  Is  discussed  at  some 
length  In  Appendix  A.  Here,  we  sviinmarlze  the  properties  which  Influence 
the  analytic  structure  of  the  eigenfunctions  and  help  to  motivate  the 
multi-component  formalism  of  Sec.  III. 

Por  paraxial  beams,  v^  Is  approximated  by 

Vr(r;Hi,Pg)  -  +  U(r;Pg)]^^^  ,  (14) 

b 


where 


U(r;Pg) 


2y^M 


-  qBAQ(r) 


(15) 
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is  the  effective  potential  that  determines  the  transverse  motion  of  a 
particle  with  anf>ular  momentum  Pq.  tfslnR  fl4),  the  perpendicular  ener^v 
can  he  expressed  as 


y  yX  +  n(r;P0). 


Oh) 


which  has  the  form  of  a  particle  Hamiltonian  in  a  time-independent 

axlsymmetrlc  potential.  Kach  particle  therefore  executes  periodic 

oscillations  in  r  while  movin.c^  aximuthally  with  an  an,(>ular  velocitv 

0  =  P./(Y,^r^).  Pxcept  in  the  special  case  where  the  net  current 
U  h 

densltv  .1  -(r)  5  J.  .(r)  +  J  -(r)  is  uniform,  the  effective  iiotential  n  is 
no  bn  pO 

anharmonic.  Because  of  this  anharmonlcl  ty ,  particles  in  <teneral  have 
precessin.^,  unclosed  orhlts  in  the  transverse  plane,  and  both  the  radial 
oscillation  frequency 


pi 

r.  '"r' 
min 


and  the  average  aximuthal  frequencv 


07) 


n,  E  <  9  > 


Vl 


r 

max 

f 

) 

’*mln 


dr 


2 

r  Iv^l 


OP) 


depend  on  H  and  P.  through  v  (r;H  ,P  )  and  the  orhlt  turnlne  points 
•L  y  r  i.  0 

^  t  P  (1*1  ’^*’6  orhlts  close  in  the  r-0  plane  onlv 

min  i  0  max  10 

when  11  /Q„  is  rational, 
r  o 

For  particles  confined  near  the  center  of  the  heam,  it  is  possible 
to  calculate  H^and  Slg  in  closed  form.  Tn  this  rei»ion,  for  an 

unhollov/ed  heam  can  he  approximated  hv 
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depends  on  both  and  P0  when  Is  not  small. 

It  Is  of  Interest  also  to  characterize  the  class  of  beam  particles 
that  have  a  particular  oscillation  frequency  since  these  particles 
are  all  resonant  with  a  wave  of  frequency  Q  >  Hg.  Properties  of  these 
particles  can  be  calculated  in  closed  form  for  the  case 

«  qSAQ(a^)  ,  (24) 

which  is  equivalent  to  the  particles  being  confined  to  r  «  a  ,  where 

n 

(19)  is  valid  and  the  effective  potential  U  is  nearly  harmonic.  We  show 
in  Appendix  A  that  v/hen  (24)  holds,  the  equilibrium  current  density 
profile  of  this  class  of  particles  is 


6Jo(r;V 


dJo(0 


26jQ(0;Qg), 


0  . 


0  <  r^<  I  R^dlg) 

I  R^dlg)  <  <  R(n0) 

R^JIq)  < 


(25) 

where  R(Jlg),  defined  as  the  maximum  radius  accessible  to  any  particle 
with  the  specified  value  of  flg,  is  given  in  the  sraall-orblt  limit  by 

R^flg)  =  3a^  (l-ingl/np)  (26) 

^2 

and  occurs  for  particles  with  |Pg|  »  4y^MR  0^/9,  which  have  orbits  midway 

between  radial  (P^  «•  0)  and  circular  (|P0|  ■  8YjjMR^JIq/9).  The  profile 

(25)  is  shown  in  Fig.  1,  and  we  note  particularly  that  5Jrt(r;no)  drops 

0  0 
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discontlnuously  to  zero  at  r  ••  ^(0^). 

Numerical  calculations  for  raonotonlcally  decreasing  profiles 

show  that  groups  of  particles  with  the  same  ftg  have  the  qualitative 
features  of  (25)  and  (26)  even  when  Is  not  small.  All  such  groups 

have  a  hollow  current  density  profile  6j^(r;nQ)  that  vanishes 
discontlnuously  at  some  outer  radius  1?(£1q).  This  radius  decreases 
monotonlcally  with  flg  because  the  pinch  force  becomes  weaker  than  that  of 
a  harmonic  potential  as  the  radius  Increases,  so  that  particles  whose 
orhlts  reach  larger  ravlll  tend  to  have  a  smaller  azimuthal  frequency  flg. 
The  plot  of  ^(^0)  In  PIR*  2  for  the  case  of  a  Bennett  profile^ ^ 

•■no''’  -  )■' 

Is  typical,  and  we  show  In  Appendix  A  that  the  curve  approaches 

P(ng)  -  1.52  (28) 

for  H,  »  qBAn(a  ).  Also,  the  particles  with  an  azimuthal  frequency  £io 
X  n  '  t7 

whose  orbits  pass  through  8(Ilg)  are  found  In  general  to  follow  orhlts 
that  are  Intermediate  between  linear  and  circular.  Since  n  and  iSlgl 
normally  have  a  one-to-one  functional  relationship,  we  can  take  Jig  to  be 
a  function  of  R  and  think  of  R  as  the  Independent  variable  characterizing 
the  class  of  particles  with  azimuthal  frequency  tlg(R).  For  later 
reference,  we  note  the  Inverses  of  (26)  and  (28): 

Jln(?)  -  +  Jlnd  -  F^/3a^),  R  «  a  (29) 

w  *“  f)  n  n 


12 


r 


n«(R)  -  +  1.52  n.a  /R,  R  »  a  . 
9  -  On’  n 


C.  Analytic  Properties  of  J(r) 

The  periodicity  of  the  particle  orbits  makes  It  possible  to  rewrite 

A 

the  Vlasov  orbit  lntej»ral  I  In  (13)  as  an  Integral  over  r  over  half  of  a 
radial  oscillation.  The  result,  derived  In  Appendix  B,  Is 


dr "‘A(r"*)  r  cos^(r)cos^(r'* ) 


I(r;H^.Pg)  -  -19c  /  (t^-)—  [ 


tan[x(nQ-  n)/n^i 


+  0(r-r')sln^(r)cos^(r')  +  0(r''-r)sln^(r'‘)co8^(r)  ]  , 


where  ^  Is  the  phase  shift  along  the  orbit, 


r  0(r';PQ)-Q 

=  /  .Ir-  -|VTP)i 


and  0  is  the  step  function. 


0(x)  -  {  ?’  ^  n 

(  I,  X  >  0. 


For  clarity,  the  and  Pg  dependences  of  the  orbit  quantities 
’^mln’  ’^inax»  ''r*  ^r’  ^0’  ^  have  been  suppressed  In  (31)  and  (32). 

A 

Although  (31)  could  be  used  In  (12)  to  calculate  J(r),  we  are  Interested 

A 

here  In  using  the  expression  to  examine  the  analytic  properties  of  J(r). 

A 

The  Integrand  of  I  is  singular  at  the  zeros  of  tan[x(n  -  fl)/n  1, 

o  r 

l.e.  for  real  JJ  and  values  of  Jlo(H,,P-)  and  0  (H,,P-)  that  satisfy 

“  X  9  r  X  o 
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0  -  0-  +  nO 
6  r 


(33) 


for  any  Integer  n.  The  principal  resonance  with  n  >  0  normally  gives  the 

strongest  coupling  because,  according  to  (32),  the  variation  of  ^  along 

an  orbit  Is  a  minimum  for  this  case,  and  the  Integrand  In  (31)  Is  least 

oscillatory.  Ue  therefore  consider  strongly  resonant  modes  with 

Re  (1  »  Im  n  and  Re  (1  -  (l.(R)  for  some  R  <  a  .  A  numerical  evaluation 

8  ~  n 

A 

of  (31)  for  realistic  A  confirms  that  the  principal  resonance  Is  Indeed 

dominant  for  these  modes.  Other  resonances,  mainly  n  -  1,  typically 

contribute  less  than  lOZ  to  the  magnitude  of  I,  while  the  nonresonant 

terms  account  for  less  than  2%.  We  are  therefore  justified  In 

substituting  ir(fla“  Q)/fl  for  tan(ir(0.-  £l)/n  ]  In  the  first  terra  of  (31) 
or  or 

and  neglecting  the  nonresonant  second  and  third  terms.  We  also  use  (23) 
to  change  variables  from  to  Wg.  These  manipulations  obviously 
preserve  the  analytic  character  of  I  neat  the  principal  resonance.  The 
perturbed  current  (12)  can  then  be  written  approximately  as 


ne(R-r) 


J(r) 


0 


dO,  I 


dP 


e 


®  iiax(0,P_)  (P^-  Pe)^^^(Pg-~PjT^^ 


X 


I  ^ 

I 


[*<')  ^  « 

N  * 


(34) 


where  the  quantities  P^(r;(2g),  given  in  Appendix  A,  are  the  maximum  and 
minimum  |Pg|  values  that  formally  give  >  U(r;Rg)  and  an  azimuthal 
frequency  +  Slg  according  to  (15)  and  (21).  The  Qg  Integration  range  In 
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I 

I 

I 


(34)  Includes  all  positive  values  of  particles  with  orbits  passing 
through  r,  and 

'^nax 

<A(r;ng,Pg)>  =  x~^n^co8^(r)  /  cos^Cr' )  (33) 

*^mln  ^ 

is  an  appropriate  orbit  average  of  A(r)  that  depends  on  flg  and  Pg 
through  r„j„.  r^^,  v^,  and 

Several  analytic  properties  of  J(r)  are  immediately  evident  from 

(34); 

(I)  In  the  special  case  "  ^(Hj-  ^^)»  where  is  a  constant, 

the  beam  current  density  profile  J^jf)(r)  is  flat  out  to  a  maximum  radius, 
and  the  pinch  potential  Is  harmonic,  provided  that  .Tpf)  is  also  flat. 

Then  flg  Is  the  same  for  all  particles,  and  when  JJ  Bg,  J(r)  becomes 
singular  as  (B  -  Bg)  ^  at  all  values  of  r. 

(II)  In  the  more  physical  case  where  ®  non-singular 

function  of  Hj^,  the  beam  profile  J^p(r)  Is  a  smooth,  rounded  function, 

A 

and  normally  dpj^(r)  Is  likewise  rounded.  Then  J(r)  has  a  singularity 
at  P(B)  for  real  values  of  B  <  B^, 

J(r)  «  In  [r  -  R(B)]  (36) 

as  r  R(B). 

A 

(III)  As  B  B^,  the  singular  point  In  J(r)  moves  to  r  ■  0.  VThen 

A 

(34)  is  rewritten  In  terms  of  It  becomes  apparent  that  J  depends  on 

r  and  B  only  through  r  =  r/R(B),  and  that  the  perturbed  field  equation 
(11)  reduces  to 


(37) 


dr  r  dr  c 

Thus  the  destabilizing  contribution  of  the  perturbed  beam  current 
vanishes  In  the  limit  Q  Aq,  and  the  hose  mode  Is  damped. 

If  we  write  (34)  In  the  form 

ng(R-r) 

J(r)  -  /  dflg  6j(r;QQ)  (38) 

to  emphasize  the  contributions  to  J  from  groups  of  particles  with  a 
common  value  of  flg,  we  see  that  5J(r;£Jg)  vanishes  at  r  “  0  and  goes 
dlscontlnuously  to  zero  at  r  =  R(Qq),  the  location  at  which 
P  (r;^fi)  ^  P  (r;^o)*  This  discontinuity  Is  responsible  for  the 
logarithmic  singularity  of  J(r)  at  r  •  R(Q). 

A 

In  the  next  section,  we  use  these  properties  of  J  to  guide  the 
formulation  of  a  macroscopic  beam  dynamics  model. 
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TIT.  Multt-Cowponent  Formalism 
A.  Bean  Dynamics  Model 

It  Is  evident  from  the  discussion  In  Bee.  TIB  that  a  beam  can  be 
thought  of  as  partitioned  Into  groups  of  particles  with  the  same  value 
of  £Jq.  Each  of  these  groups  Is  localized  within  <  r  <  B(ng),  and  all 
particles  In  such  a  group  are  resonant  with  a  wave  of  frequency  ft  ■  ftg. 

A 

This  radial  localization  leads  to  the  slngiilarlty  In  J(r)  at  R(ft0)  found 
In  Bee.  TIC.  A  macroscopic  bean  dynamics  model  might  therefore  represent 
the  resonant  structure  correctly,  provided  that  It  localizes  the  part  of 
the  bean  resonant  at  frequency  ft  within  radius  B(ft). 

Our  approach  Is  to  represent  the  beam  slice  at  axial  location  (  as  a 
superposition  of  components  with  different  radii,  each  with  an 
equilibrium  current  density  profile  of  the  form 

i6j-(0;R)0(r/l>),  r  <  R 

0,  r  >  R, 

where  0(0)  «  l.  The  quantity  R  Is  used  here  to  identify  the  component 
with  radius  R,  just  as  either  ftg  or  R  was  used  In  Sec.  IT  to  Identify  the 
group  of  particles  with  the  same  ftg.  The  bean  current  d|jQ(r)  is  the  sum 
of  the  component  densities, 

Q» 

Jjj^(r)  -  /  dR  6J^(r;R).  (40) 

To  treat  the  linearized  beam  response  to  an  ra  ■  I  wave  of  the  form 

A 

(10),  we  specify  that  J(r)  be  the  sun  of  the  Individual  component 
perturbed  densities. 
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(41) 


A  A 

J(t)  -  /  dR  6j(r;R), 

0 

and  we  require  that  each  component  respond  rigidly  to  the 
linearized  J  x  R  force  averaged  over  Its  cross  section,  so  that 


«J(r;R) 


-  C(R)  6jQ(r;R) 


(42) 


where  C(R)  Is  the  amplitude  of  the  Fourier-analyzed  transverse 
displacement.  From  Newton's  force  law,  we  find  after  two  Integrations  by 
parts  that 


-n  C(R) 


xqB 

Yj^M«I(R) 


/  ‘Ir  r[«J^(r;R)(^  + 


dA. 

+  6j(r;R)  -jj— ], 


(43) 


where 


R 

6I(R)  -  2it|  dr  r«Jrt(r;R) 


(44) 


Is  the  component  current.  We  note  that  the  second  bracketed  term  In  (43) 

represents  the  restoring  force  on  the  displaced  component  due  to  the 

equilibrium  field.  Using  (42),  this  term  can  be  written 

as  -fl^(R)C(R),  where 
c 


n^(R) 

C 


xqg 

Tj^M6I(R) 


F 

/  dr  r 

n 


**^0  3 

jj-  5J„(rilO 


(45) 


Is  the  natural  oscillation  frequency  of  the  component  with  radius  R  in 
the  equilibrium  field.  Kquatlon  (43)  can  then  be  rewritten 


i(F) 


xqg 

Y,jM«I(R) 


1 

Sl^(R)  - 
c 


/  dr  rA(r)  6jj^(r;R), 
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(46) 


ilii 


y-  ■  ■  " 


It  Is  evident 


2  2 

which  explicitly  displays  the  resonance  at  H  (R)  ■  n  . 

c 

that  plays  a  role  similar  to  that  of  the  individual  particle 

frequency  i2g(R),  although  we  should  recall  that  the  component  executes 
linearly  polarized  simple  harmonic  oscillations  with  no  azimuthal  motion, 
whereas  Individual  particles  follow  complicated  precesslng  orbits  in  the 
transverse  plane. 

The  Vlasov  analysis  of  Sec.  II  shows  that  the  group  of  particles 
with  a  given  value  of  does  not  in  fact  respond  rigidly  to  a  transverse 
force.  Nonetheless,  we  can  choose  the  component  shape  (?(r/R)  so  that  the 

A 

perturbed  current  density  6j(r;R)  of  a  component  reproduces  the 

A 

properties  of  the  Vlasov  perturbed  current  contribution  6J(r;ftQ).  As 

A 

noted  in  Sec.  IIC,  6j(r;R)  should  be  acontinuous  function  on 
n  <  r  <  R,  that  vanishes  smoothly  at  r  =  0  and  drops  dlscontlnuously  to 
zero  at  r  =  R.  We  choose  the  convenient  mathematical  form 

6j^(r;R)  =  6Jf^(n;R)  (1  -  rV^)  OCR  -  r),  (47) 


which  satisfies  these  requirements  and  yet  allows  the  component  profile 
to  be  varied  from  a  triangular  shape  to  a  step  function  by  varying  A  from 
unity  to  infinity. 

With  this  component  shape,  (40)  can  be  solved  for  the  coefficient  of 


«Jf,(0;R) 


7 

d'J 


bO 


dr 


dJ 


(1-1) 


bO 


dr 


Ir- 


(4«) 


and  the  component  current  (44)  is 
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(49) 


«I(R)  -  ^  . 


Combining  these  expressions  with  (41).  (42),  and  (46)  then  gives 


dX 


,  (£+2)qe  t-1  ,  dR  r_  bn  ”''bO  , 

^  J  —T  -  -dF-  JrH 


r  R 
1 


dr 


/  dr'r'^  A(r'), 


fl^(R)  -  n 
c 


(50) 


where,  from  (45), 


fl.(R) 


4ffq6  <d-2  1 


V'  ‘  „ 


(51) 


Comparing  (51)  with  (22),  we  see  that 


«'(">  -  «0  * 


(52) 


Independent  of  t.  For  real  JJ  values  less  than  the  maximum  0^,  the 
resonant  denominator  In  (50)  can  be  represented  near  the  singularity  as 


dn"' 

"  -  'cW  (  # 


^R-R  (11)  • 
c 


where  F  (il)  Is  the  radius  of  the  component  with  fl  ■  Jl  and  may  be  thought 
c  c 

of  as  the  Inverse  function  of  (51).  This  pole  In  the  Integrand  of  (50) 

leads  to  a  logarithmic  slngxilarlty  In  J(r)  at  r  -  P  ((1),  In  agreement 

c 

with  the  Vlasov  result. 

To  choose  an  appropriate  value  for  the  Index  t  In  the  component 

A 

profile  (47),  we  examine  the  .T(r)  expression  (50)  as  0  ♦  Jl^.  Since 
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o 

n  -  il  (R)  varies  like  R  for  small  R,  the  resonant  denominator  In  (50) 
0  c 

can  he  approximated  for  fl  Jl.  by 


We  can  then  write  (5^^)  in  terms  of  r  =  r/R  (0)  and  R  =  R/P  (fi): 

c  c 

J(r)  -  -  ^  I- -  /  — ^  --  -  f  dr'r'^A(r').  (53) 

V’  R\n)  -  R^-i  fi 

(I  c  r 

We  note  from  (4ft)  that  for  Z  =  2,  6Jp|(r';R)  “  P  for  sp«11  R,  whereas  for 

A 

any  other  value  of  Z,  6J|^(0;R)  «  r.  Consequently,  as  P^(0)  ■*■  0,  J(r)  Is 

Independent  of  R  (0)  when  Z  =  2  and  becomes  negligible  compared  with 
c 


dr  r  dr 


R  (0)  dr  r  dr 
c 


so  that  the  A  equation  (11)  takes  the  form  of  the  Vlasov  result  (37). 

''  -2 

For  other  values  of  4,  J  *  R  (W)  in  the  same  limit.  We  thus  conclude 

c 

that  A  =  2  Is  the  only  physically  acceptable  choice. 

With  £  =  2,  the  multi-component  expression  (53)  for  6J  also  shows 

remarkable  agreement  v;lth  the  Vlasov  result  In  other  respects.  The 

multi-component  resonance  location  P^(Sl)  Is  compared  with  the  Vlasov 

resonance  location  R(£J)  in  Fig.  2.  The  two  quantities  agree  to  within  5% 

everywhere  and  to  much  greater  accuracy  for  the  large  fraction  of  beam 

particles  with  r  /a  <  2.  In  the  limit  JJ  (51)  gives 

max  n  ~  n’ 


a^(P)  -  a„(i  - 
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for  A  ■  2,  in  exact  agreement  with  the  Vlasov  expression  (29)  for  n^CR). 
In  the  opposite  limit,  0,  the  multi-component  expression  is 

n^(R)  -  2^"^^  Vn^*^* 

the  coefficient  here  agreeing  with  the  Vlasov  result  (30)  to  within  5%. 

In  addition,  a  numerical  evaluation  of  J(r)  for  the  two  models,  using 
realistic  J^jQ  and  A  profiles  in  (12)  and  (50),  shows  agreement  within  7% 
in  the  perturbed  current  magnitude.  This  suggests  that  the  multi- 
component  model  accurately  represents  the  strength  of  resonances,  as  well 
as  their  location  and  analytic  form. 


B.  Calculation  of  Elgenmodes 

We  formulate  a  hose  eigenvalue  equation  by  substituting  the  multi- 

A 

component  J  expression  (50)  with  =  2  into  the  m  »  2  field  equation 
(11).  The  result  is 


did  : 

-z - rA 

dr  r  dr 


4x10)0 

2 

c 


A 


00  H  T 

16xqB  r  ^  (  L.  1  pO  )  1 _ 

r  r"^  n^R)  -  0 

C 


/  dr'r'^  A(r'), 
0 


(57) 


where,  from  (51), 


(58) 


Boundary  conditions  are  obtained  from  the  requirement  that 
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A(0)  -  0 


(59) 


for  an  m  °*  1  mode  and  from  the  far-field  solution  of  (57) 


A 


(60) 


This  eigenvalue  problem  can  be  solved  for  either  (>>(0)  or  0((i)).  It  Is 
apparent  from  (10)  that  a  wave  with  real  0  and  Im  u  >  0  has  constant 
amplitude  at  a  given  beam  slice  but  grows  spatially  back  from  the  beam 
head  (i.e.  with  Increasing  c).  Similarly,  a  mode  with  real 
u  grows  with  Increasing  z  If  Im  0  >  0.  We  concentrate  here  on  modes  with 
real  0  because  these  modes  show  the  effects  of  localized  resonances  most 
clearly. 

To  solve  the  dispersion  equation  numerically  we  first  rewrite  (57) 
as  a  homogeneous,  fourth-order  ordinary  differential  equation, 


d  I  d  1  d^  /  dr  -_•> 
dr  r  dr  r  dr  dr  ^ 


Airlut 


dr 


2r)] 


16nq0  _J.  r  ^  1  ^^60  ^  F 

Y  Me  ^2  ^  dr  r  dr  ^^2^^^  _  „2 

c 


(61) 


where 


r(r)  S  4  /  dr'r'^  A(r').  (62) 

r  0 

A 

This  definition  of  T  and  the  two  constraints  on  A,  (59)  and  (60),  give  an 
appropriate  set  of  boundary  conditions: 
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(63) 


r(o)  -  0 

r'(n)  -  0 

115  «  r'(r)  -  0 

li^i  ^  r"(r)  »•  0. 

We  reduce  (61)  to  a  homogeneous  matrix  equation  by  representing;  the 
derivatives  by  finite  differences  on  a  nonlinear  mesh,  with  prld  points 
clustered  near  the  axis  and  any  resonances.  For  real  appropriate  .lump 
conditions  at  resonance  points  are  provided  by  local  power  series 
solutions  of  (61).  Kl{;enf requencies  are  calciilated  by 

specifying  a>  or  fl  and  iteratively  ad  .justing  the  other  frequency  until  the 

matrix  determinant  vanishes.  The  corresponding  eigenfunctions  are 

1  ? 

calculated  by  the  method  of  Inverse  Interation. 


C.  Equilibrium  Specification  for  Beams  with  Beturn  Current 

Although  the  dispersion  equation  (61)  can  be  solved  for  any  choice 

2 

of  JpQ(r)  and  Jjjf^(r)  that  give  a  non-negative  an  arbitrarily  chosen 
may  correspond  to  a  physically  unlikely  distribution  function.  Here, 
Instead,  we  specify  a  Hax\7ell-Boltzmann  distribution  function. 


"  2Tr'Y;Mq8cHj^  exp(-  —  ), 


(6A) 


and  also  specify  the  conductivity  profile  a(r)  and  the  effective  net 
current  fraction,  defined  as 


f 

e 


2xc 


/  dr  r^U 
0 


dA 

bn 


0 


(65) 
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The  equilibrium  electric  field  E^q  *  -SdA^/dC  Is  taken  to  be  Independent 
of  r,  which  Is  a  good  approximation  for  r  <  a>  The  equilibrium  beam 
current  expression  (9)  can  be  evaluated  analytically  In  this  case  and 
gives 


°  exp(qBAQ/H^^), 


(66) 


and  substituting  this  expression  Into  (65)  gives 


f 

e 


2c  „ 
q6Ib  »b‘ 


(67) 


an  expression  of  the  Bennett  pinch  condltlon^^.  For  a  given  the 
effective  net  current  f^  depends  only  on  the  perpendicular  temperature 
and  can  be  specified  Independently  of  JpQ.  The  equilibrium  potential  Aq 
Is  calculated  by  numerically  Integrating  Ampere's  equation  (2)  In  the 
form 


1  d  r  T  /  \  /  w  1 


(68) 


where  JjjQ  Is  given  by  (66)  and  (67)  as 


e  b 


(69) 


and  E^q  Is  determined  Iteratively  by  requiring  that 


-  2n  /  dr  r  J^^Q(r). 


(70) 
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This  procedure  yields  a  Bennett  current  profile 


t 

11 


-2 


(71) 


with  a  beam  radius 


>1/2 


(72) 


either  when  E^q  ~  0  or  when  a  has  the  same  radial  dependence  as  (71 )• 


D.  Comment  on  the  "Energy-Group"  Model 

In  Ref.  6  we  used  a  multi-component  model  with  a  simple  component 
current  density 


djQ(r5R)  -  djQ(0;R)  0(R-r), 


(73) 


which  corresponds  to  (47)  with  1  These  components  were  referred  to 

as  energy  groups  because  each  one  represented  the  group  of  particles  with 
some  particular  value  of  perpendicular  energy  Unfortunately,  the 

present  analysis  shows  that  this  component  shape  Is  Incorrect.  The 
result  of  using  (73)  In  the  multi-component  J  expression  (50)  Is 


J(r) 


_gB_ 

V 


A(r) 


n"(r)  - 

c 


dJ 


bO 


dr 


(74) 


Equation  (74)  gives  a  pole  In  J(r)  at  r  ■  R^(n),  rather  than  the 
logarithmic  singularity  found  by  the  Vlasov  analysis  of  Sec.  II. 
Substituting  (74)  Into  (11)  leads  to  a  second-order  dispersion  relation 


r 


^  1  rA(r)  +  A(r) 

dr  r  dr  2 

c 


Airgg  A(r)  j_  **‘^bO 

n^(r)-  ^ 


(75) 


which  duplicates  (45)  in  Pef.  6,  In  the  limit  (2  ♦  Q^,  (75)  can  be  put  in 
the  form 


.  ,  ,  «  4iria)P^  (£2)0(0) 

—  -  —  rA(r)  + - ^ - 

dr  r  dr  c^ 


A(r) 


■'bo'"'  * 

^o'"'  I 


2 

n 

2  ’ 


(76) 


where  r  =  r/R  (£2),  a_  is  the  net  current  scale  radius  defined  hy  (15), 

C  u 

and  a  is  the  analogous  beam  current  radius.  According  to  (76), 

2  —2  -1 
uR  (£2)  remains  constant  as  £2  +  £2*.  Rince  R  (£2)  «  ($2^“  ®)  >  tills  model 

c  »l  c  t* 

predicts  a  pole  in  u)  for  £2  *  £2^  and  an  infinite  growth  rate  Im  u 

as  £2  ♦  £2^  from  below.  The  correct  result  is  that  the  wave  is  damped  in 

this  limit. 

This  sensitivity  of  the  solution  to  the  component  shape  motivated 
the  Vlasov  analysis  of  5ec.  IT  and  the  more  careful  choice  of  component 
shape  In  Sec.  TIIA. 
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IV.  Results 


We  have  numerically  solved  the  multi-component  hose  eigenvalue 
problem  discussed  In  Sec.  Ill  for  real  ^  and  various  choices  of  effective 
net  current  fraction  f^,  beam  current  profile  conductivity 

profile  a(r).  The  results  are  discussed  In  this  section. 

A.  Eigenfunctions 

(I)  Resonant  eigenfunctions.  In  the  usual  case  where  the  net 

current  density  profile  bell-shaped  and  peaked  on  the  axis, 

n  (R)  also  peaks  on-axls  and  modes  with  real  <  0.  are  resonant  at 
c  u 

r  >  R^($l).  A  typical  hose  eigenfunction  for  such  a  mode  Is  shown  In 

Fig.  3.  As  expected,  there  Is  a  logarithmic  singularity  In  the  perturbed 

current  at  r  >  R  (ii),  and  we  note  that  the  perturbation  Is  essentially 
c 

confined  within  0  <  r  <  R^.  Varying  Q  or  f^  changes  R^.  and  therefore  the 
width  of  these  resonant  eigenfunctions,  but  their  shape  Is  hardly 

A  A 

affected:  If  A  and  J  are  plotted  as  functions  of  r/R^.,  they  are  nearly 
Independent  of  Jl  and  f^  so  long  as  Is  monotone  decreasing  with  R. 

The  resonant  structure  of  these  elgenmodes,  properly  regarded  as 
Internal  disruptions  localized  within  resonances.  Is  Ignored  In  the 
earlier  rigid-beam^ and  spread-mass^  hose  dynamics  models,  where  the 

A 

perturbed  current  J(r)  Is  taken  to  be  simply  a  rigid  displacement  of  the 
entire  beam. 

(II)  Doubly  resonant  eigenfunctions.  When  the  radial  profile  of 

JpQ(r)  Is  narrower  than  that  of  large  enough,  the  net 

current  density  peak  away  from  the  axis.  In  this 

case,  (2^(R)  also  peaks  off-axis,  and  a  mode  with  real  U  <  max(n^)  Is 
resonant  at  two  different  radii  that  satisfy  ^^(R)  “  We  then  find 
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Fig.  3  —  (a)  Real  and  imaginary  parts  of  A(r)  for  a  typical  hose 
eigenfunction  with  n  real.  In  this  case,  J|,q  and  o  have  Bennett 
profiles  with  radius  a,  and  we  have  n/flg  =  0.8  and  (f2)/a 

»  0.7.  (b)  Real  and  imaginary  parts  of  J(r)  for  the  same  mode 
as  (a). 
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that  the  perturbation  Is  largely  confined  between  the  two  resonant  radii, 
as  shown  In  Fig.  4. 

(Ill)  Mon-resonant  eigenfunctions.  Modes  with  or  with 

significant  In  are  not  resonant.  We  find  then  that  the  form  of  the 
eigenfunctions  Is  closer  to  that  of  a  simple  displacement  of  the  beam, 
where 

J(r)  «  dJ^^/dr. 

(iv)  Higher  radial  elgennodes.  ’’'or  any  value  of  ft,  there  Is  an 
Infinite  sequence  of  elgennodes  corresponding  to  Increasing  radial 
quantum  number.  ’’’or  beams  with  a  Bennett  profile  and  no  return  current, 
we  Invariably  have  found  these  modes  to  be  strongly  damped.  Farller  work 
Indicates  that  some  of  these  nodes  may  be  unstable  for  flat  beam 
profiles  and  for  cases  with  a  large  return  current  fraction  »  ,  but  we 

have  not  examined  these  cases. 


B.  Plsperslon  Characteristics 

The  eigenvalue  equation  (57)  Is  conveniently  put  into  dimensionless 
form  by  Introducing  the  scaled  quantities  r/a,  »  and  ft/ft^l^,  where 

T  (77) 

Is  the  decay  time  for  dipole  flelds^’^  and 


2lTq6Jj^^(0) 


(78) 


2 

Is  the  value  ft^^  would  take  In  the  absence  of  return  current,  given  In 

1  2  2 

terms  of  the  beam  plasma  frequency  u  (r)  by  -r-  B  w  (0).  The  resulting 

p  £  p 
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eigenvalue  problem  is  Independent  of  a,  6,  y^,  o(0),  and  The 

dispersion  relation  can  therefore  be  quite  generally  expressed  as  a 

functional  relationship  between  and  depending  only  on  the 

scaled  quantities  a(r)/o(0),  J^^Cr )/Jjjfj(0) ,  and  present 

our  dispersion  results  In  this  form. 

(1)  No  return  current.  The  hose  growth  rate  J.m(u)T,)  and  the 

d 

frequency  Re(u)T.)  are  plotted  as  the  solid  curves  In  Fig.  5  for  the  case 
d 

where  J^jf^(r)  and  a(r)  have  Bennett  profiles  of  equal  width.  We  note  that 

Instability  occurs  only  over  the  range  0  <  ^  0.7B,  with  the  maximum 

growth  rate  ojt,  ••  0.65  occurring  at  -  0.57.  XJe  note 

d  DO 

that  u)  ■  0  when  J)  •*  0,  so  that  a  beam  displaced  rigidly  from  the 
conductivity  channel  axis  la  neutrally  stable,  reflecting  the  fact  that 
there  is  no  net  force  on  the  beam  in  the  absence  of  plasma  current. 

For  0  <  <  1  the  growth  rate  is  determined  by  competition  between  the 

resonant  wave-particle  coupling  and  phase  mixing  due  to  the  betatron 
frequency  spread.  Fven  though  resonance  effects  are  stronger 
as  fl  the  number  of  resonant  particles,  given  approximately  by  (49), 

decreases  as  -►  0,  and  for  >  0.78  the  damping  dominates.  A 

dispersion  relation  of  this  type  leads  to  convective  instability  In  the 
beam  frame^,  l.e.  a  perturbation  introduced  at  some  point  along  the  beam 
will  grow  as  It  propagates  back  In  C  but  will  eventually  decay  at  any 
given  value  of  C. 

This  case  has  previously  been  studied  using  the  energy-group  raodel^ 
and  the  spread-mass  model An  exact  solution  of  the  energy-group 
dispersion  equation  (75)  gives  a  growth  rate  Tm  u(n)  that  agrees  well 
with  our  present  result  for  ^/^^o  ^  ^ut  continues  to  Increase  to  a 

pole  at  1,  for  reasons  discussed  In  Bee.  Tlin.  However,  the 
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variational  method  used  In  Ref.  6  to  approximately  solve  (75)  reduces  the 

resonance  strength  and  leads  to  the  dashed  curves  In  FIr.  5,  which  agree 

well  with  the  multi-component  result.  In  the  spread-mass  model,  each 

slice  of  the  beam  Is  subdivided  into  rigid  components  which  have  density 

profiles  Identical  to  the  beam  Itself.  The  perturbed  current  Is  thus 

constrained  to  be  a  rigid  displacement  of  the  slice  as  a  whole,  and  there 

Is  no  correlation  between  the  radial  extent  of  a  component  and  Its 

resonant  frequency.  The  components  carry  eaual  fractions  of  the  beam 

current,  but  a  realistic  spread  in  their  oscillation  frequencies  Is 

Introduced  by  assigning  the  components  different  masses.  For  a  beam  with 

a  Bennett  density  profile,  the  resulting  dispersion  relation  Is  Identical 

to  the  variational  approximation  of  the  energv-grotin  expression,  plotted 

In  Fig.  5  as  dashed  curves,  and  It  agrees  remarkably  well  with  our  multi- 

component  result  over  the  entire  unstable  range.  "Hie  reason  for  this 

good  agreement  Is  that,  for  a  Bennett  profile  with  no  return  current,  the 

unstable  frequencies  have  resonant  radii  R  (il)/a  >  1,  so  that  the 

c 

components  which  Interact  most  stronglv  with  the  wave  are  nearly  as  broad 
as  the  spread-mass  components.  Consequently,  the  resonant  components  In 
the  two  models  respond  similarly.  For  waves  with  ~  1 

and  R/a^«  1,  which  have  stronger  resonant  effects,  the  two  growth 
spectra  become  quite  different,  even  though  the  modes  are  damped  and 
therefore  physically  unimportant. 

The  growth  spectra  for  both  the  spread-mass  and  multi-component 
models  are  sensitive  to  the  equilibrium  beam  density  profile  •^bn  ,  and  the 
two  models  do  not  agree  as  well  for  other  beam  profiles  as  they  do  for 
the  Bennett  profile.  For  a  Causslan  profile 
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Fig.  5  —  Real  eigenCrequencies  Re  u>  and  growth  rates  Im 
u  for  the  hose  inodes  with  n  real,  calculated  for  a  Bennett  beam. 
The  solid  curves  are  obtained  using  the  multi-component  model; 
dashed  curves  come  from  the  spread-mass  model  or  the  variational 
approximation  to  the  energy-group  model. 
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(79) 


J^^(r)  “  exp(-r^/a^) 


the  spread-mass  formalism  can  be  applied  onlv  If  the  conductivity  has  the 
broadened  form 


a(r)  =  a(n)  ^ 


■’  2 

e.xp(-r  7a  ' ) 


1  -  exp(-r  /a  ) 


The  growth  spectrum  from  the  multi-component  treatment  Is  shown  for  this 

case  as  the  solid  curve  in  Fig.  b.  The  curve  Is  more  sharplv  peaked  near 

the  maximum  growth  rate  than  the  Bennett  result  in  Fig,  S  because  there 

are  fewer  particles  with  near  0  and  1.  The  corresponding  spread- 

c  nu 

mass  result,  plotted  In  Fig.  6  as  a  dashed  curve,  is  broader  than  the 
mtiltl-conponent  curve  and  Is  shifted  to  larger  SI  values  since  the  mass 
distribution  In  this  case  has  more  low-mass  particles.  The  dotted  curve 
In  Pig.  6  is  the  growth  spectrttm  calculated  from  the  mviltl-component 
dispersion  equation  using  (79)  for  and  a  similar  Gaussian  profile 
for  0.  This  is  a  more  realistic  conductivity  profile  than  (FO)  for  a 
Gaussian  beam  and  gives  a  substantially  higher  maximum  growth  rate 
beca\ise  the  wave-driven  dlnole  return  current  is  confined  nearer  the 
axis.  This  sensitivity  to  the  width  of  the  conductlytty  channel  Is  a 
general  feature  of  the  hose  instability  and  is  dlsciissed  subsequently. 

(11)  Effect  of  return  current.  The  presence  of  eouillbri>im  plasma 
current  worsens  the  Instability  hy  magnetically  repelling  the  displaced 
beam,  and  this  repulsion  Is  Increased  when  a  narrow  o  channel  confines 
JpQ  and  the  perturbed  plasma  current  near  the  axis.  Plg\ire  7  shows  the 
effect  on  the  hose  growth  spectrum  of  a  nonzero  return  current  fraction 

f„  =  1-f  when  J,„  and  a  both  have  Bennett  profiles  with  the  same 
r  e  bn 
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Fig.  6  —  Hose  growth  rates  for  a  beam  with  a  Gaussian  current 
density  profUe.  The  multi-component  model  gives  the  dotted 
curve  for  a  Gaussian  conductivity  proflle  and  the  solid  curve  for 
a  broadened  a  profile.  The  dashed  curve  comes  from  the  spread- 
mass  model,  which  required  the  broadened  a  profile. 


radius.  The  return  current  Is  seen  to  be  the  dominant  destablllzlnR 


factor  when  >  0.3,  and  for  these  f^  values  the  maxlnum  growth  rate 
Increases  approximately  proportionally  to  f^.,  as  expected  for  pure 
magnetic  repulsion.*’  In  this  case,  Im  o)  >  0  even  for  0  =  0,  so  that  a 
beam  that  Is  rigidly  displaced  from  the  conductivity  channel  axis  Is 
pushed  further  away.  This  Is  an  absolute  Instability  In  the  beam  frame. 

(Ill)  Narrow  conductivity  channel.  A  narrow  conductivity  channel 
can  also  substantially  Increase  hose  growth  by  preventing  the  perturbed 
plasma  current  from  following  the  beam  distortion,  even  when  there  Is  no 
equilibrium  plasma  current.  The  growth  spectra  in  Tig.  R  are  calculated 
for  beams  with  V"  and  a  bennett  profile  of  o(r)  with  a  radius 
different  from  that  of  the  beam.  We  see  a  dramatic  Increase  In  the 
maxlmttm  growth  rate  when  the  ratio  a  of  the  conductivity  radius  to  the 
J^Q  radius  Is  less  than  unity.  Broadening  the  conductivity  channel,  so 
that  a  >  1,  has  a  stabilizing  effect  on  hose  modes,  but  the  effect  Is 
less  pronounced  than  for  narrow  channels  because  A  is  still  confined 
within  and  largely  determines  the  radial  profile  of  the  plasma 

current  when  a  is  broad  and  =  0. 

(iv)  Return  current  and  narrow  conductivity  channel.  Figure  9  shows 
the  growth  spectra  for  a  case  vrlth  equilibrium  return  current  f^.  “  0.23, 
and  several  a  values.  As  expected,  the  presence  of  return  current  and  a 
narrow  conductivity  channel  act  together  to  Increase  the  growth  rate  even 
more  rapidly. 

C.  Concluding  Remarks 

We  have  developed  a  tractable  model  of  the  resistive  hose 
Instability  which  retains  many  features  of  a  full  Vlasov  treatment,  and 


we  have  used  the  model  to  examine  the  dispersion  relation  and 
eigenfunctions  in  a  variety  of  cases.  The  rigid  component  dynamics  used 
here  is  particularly  suited  to  m  ••  1  perturbations  like  the  hose 
instability,  but  with  appropriate  modifications  this  approach  is 
applicable  to  m  =  n  beam  modes  and  possibly  to  other  resonant 
instabilities  of  periodic  systems. 

We  note  that  these  hose  stability  results  are  restricted  to  beams 
propagating  in  a  medium  with  a  fixed  conductivity  profile,  i^or  the  case 
of  Intense  beams  propagating  into  initially  un-ionized  or  weakly  ionized 
gas,  the  ionization  and  heating  generated  by  the  beam  typically  causes 
the  conductivity  channel  to  follow  beam  displacements.  Since  this  effect 
reduces  the  separation  between  the  beam  and  plasma  currents,  it  is 
stabilizing.  We  have  shown  in  the  related  case  of  the  sausage 

instability^  that  this  effect  can  be  treated  correctly  only  by  including 
conductivity  generation  consistently  in  the  equiltbriun  as  well  as  in  the 
perturbation,  l.e.  the  equllihrlum  conductivity  a(r,C)  must  he  allowed  to 
increase  with  This  self-consistent  treatment  of  beam-generated 
conductivity  Introduces  non-ignorable  C  dependence  and  Is  thus  beyond  the 
scope  of  the  one-dimensional  eigenvalue  treatment  developed  in  this 

O 

paper.  However,  we  have  developed  a  numerical  treatment  in  which  the 
multi-component  model  is  solved  by  following  the  evolution  of  an  initial 
perturbation  in  space  and  time,  rather  than  by  Fourler-analyztng  in 
C  and  z.  This  approach  allows  us  both  to  treat  conductivity  generation 
self-consistently  and  to  Include  axial  variations  in  the  enullihrlum  beam 
and  return  currents.  Since  we  do  not  assume  high  conductivity  or  charge 
neutrality  in  this  time-dependent  model,  we  can  also  treat  the  weakly- 
pinched  beam  head.  This  work  will  be  reported  in  a  separate  pvibllcatlon. 
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It  Is  also  to  be  noted  that  tn  cases  of  very  high  beam  current 
density  and  low  enough  background  gas  density,  such  as  tn  ton  beam 
transport  for  Inertial  confinement  fusion,  heating  of  the  background  gas 
can  lead  to  significant  hydrodynamic  notion  of  the  gas.^^  This  effect 
should  also  be  Included  self-conslstently  in  a  full  treatment  of  hose 
instability  In  such  systems. 
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Appen<!lx  A;  Particle  Orbit  Properties 
(1)  Particle  motion  In  the  transverse  plane 

Tlsings  the  definitions  (3)-(5),  the  energy  R  of  a  particle  In  a 
high-conductlvlty  beam  can  he  written  In  terras  of  the  radial  momentum 
component  p^  and  the  canonical  momentum  components  Rg  and  Pj,  as 

H  =  yMc^-  +  c^[  Pr+-f  + 


Equation  (Al)  may  be  regarded  as  an  equation  for  p^  as  a  function  of  r 

and  the  three  constants  of  motion,  ^or  a  paraxial  beam 
2  2  2 

with  p  »  p  +  Po  and  v 


Ubl  «  =  (Y^-l)^^^Mc^|q^^ 


(A2) 


and  from  Ampere's  equation  (2)  we  find 


IqA^/cI  <  Iql^/c^l  «  lol^/c^l  •>  YMc  -  p^. 


(A3) 


This  allows  us  to  rewrite  (Al)  In  the  approximate  form 


r 


(A4) 


where 


2  ^  „2  A 
(Pc  +  M  c 

'■  z 


1/2 


(A5) 


Is  the  energy  associated  with  P^. 


Since 


2 


Is  large  compared  with  the 
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energy  In  perpendicular  motion  and  In  (A4)  can  be  approximated  by 


7  1  2  ^0  '’’’z 

H  -  y^Mc  +  (p^  +  -5-)  -  A^. 

D  r  D 


(A6) 


The  radial  velocity  v  =  p  /(yfO  “  P  /(Y.W)  can  then  be  written 
■  r  r  r  D 


v2(r;«i.Pe)  [h^  -  iKrjPg)]. 


(A7) 


where 


=  H  -  y^nc\ 


(A8) 


may  be  regarded  as  the  energy  In  the  transverse  plane,  and  IT  Is  an 
effective  potential  given  bv 

1  ^0 

n(r;P0)  s  _-q3A^(r).  (AP) 

'b‘  r 

The  complete  particle  dynamics  thus  reduces  to  motion  In  the  transverse 
plane  subject  to  an  effective  central  potential. 


(11)  Particle  oscillation  frequencies 
2 

Since  v^  Is  a  function  of  r  only,  for  a  given  and  Pg,  particles 

execute  periodic  radial  oscillations  while  gyrating  at  an  angular 

•  2 

velocity  6  ••  ^0/TjjIlt  •  The  characteristic  frequencies  are  the  radial 
oscillation  frequency 

r 

max  ,  , 

V«i.p0)  - (Ain) 

'‘min  ^ 
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and  the  mean  azimuthal  frequency 


0  (H|,P£,)I*g  *^nax  • 


(Ain 


r  .  r  |v  I 
min  r 


where  the  turning  points  r  .  (Hi.Pa)  *■  **■*  solutions  of 

iQxn  jl  V  nftx  X  o 


v^(r;H^.P0)  -  0 


(A12) 


such  that  H,  >  U(r;P«)  for  r.<r<r  .  In  general,  ft  /ft.  Is 
J.  min  max  r  u 

Irrational,  and  the  orbits  do  not  close.  In  the  special  case  of  a 

harmonic  potential,  corresponding  to  a  uniform  net  current,  these 

frequencies  are  Independent  of  H,  and  P.,  and  ft  -  2ft.. 

X  0  r  0 

For  particles  confined  near  the  axis,  where  the  potential  Is  nearly 
harmonic,  we  can  express  ft^,  and  ftg  In  closed  form.  In  this 

limit,  the  net  current  can  be  expanded  In  the  form 

•'no"'  •  <«’> 


and  A^Cr)  can  be  written  to  the  same  accuracy  as 

,  7  2 

ArjCr)  -  -  f  (1  - 

2a 

n 

where  a^  Is  the  radial  scale  length  of  the  net  current  profile  and  we 
assume 


«  1. 


(A15a) 
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which  holds  for  particles  with 


«  qflAQ(a^)  . 

Using  (A14),  we  can  reduce  (AlO)  and  (All)  to 


■  - 


+  2r^ 

max  min 


2a 


2  2 
r  “  r  . 
max  min 


2a^-  r^  -  2r^ 
n  max  min 


)] 


-1 


and 


a  P. 
r  n  0 


0  xQ, 


21/2 


r%  [2a^-  r^  -  2r%  ) 
mln^  n  max  min-' 


1/2 


n(- 


2  2 
^mln  *^max 


min 


r^  -r' 
_ max  min 

2a^-  r^  -  2r^ 
n  max  min 


), 


where  K(o)  and  n(a,0)  are  complete  elliptic  Integrals  of  the 
third  kinds 


K(a) 


1 

^  I 
0 


dx 


[(1-X^)(1-(K^ 


)] 


1/2 


a  «  1 


and 


(A15b) 

(A16) 

(A17) 

first  and 

(A18) 
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n(a,0)  =  /  j—  ^  2  1 1/2 

- S72  -  'Ifs  1  • 

(1-$)  ' 


(Al<)) 


a  «  1  . 


In  (A16)  and  (A17)  we  have  Introduced  the  characterlatic  frequency 


n  =  r  2131  j  (0)1^^^ 

0  ■  ^  YwMc  * 

D 

and  a$;ain  r  .  and  r  (P|»Pq)  are  defined  by  (A12). 

rain  J.  o  max  x  0 

order  In  r^^/a„. 


To  zeroth 


9  9 


v4 


(A2n) 


Using  this  approximation  and  the  small  argument  expansions  of  K  and  11 
in  (AlA)  and  .(A19),  the  orbit  frequencies  (A16)  and  (A17)  can  be 
rewritten  as 


’r  ■  -  —2:1  ) 


(A21) 


3U,- 

«0  -  8gn(Pg)n„(l  -  2  2  -  ). 


We  see  that  Is  the  magnitude  of  In  the  limit  0, 


(A22) 


Tn  the  opposite  limit 


(A23a) 


»  1, 

max  n 


which  is  equivalent  to 


(A?3b> 


we  can  make  the  approximation 


A-(r)  *  Ac  Sl„l  ln(r/a  ) 
u  n  n 


over  most  of  a  particle  orbit.  Writing  In  terras  of  the  outer  turning 

point  r  (H.,0)  of  a  particle  with  P.=  0,  we  then  have 
max  1  D 

.2 

2^2 


"i-  “"’"s'  ■  V'«^^  - 

b  r 


(A2A) 


and  the  orbit  frequency  expressions  (AlO)  and  (All)  reduce  to 


fl  - 


a  Qf. 
n  0 


r  -1/2  r  (H,,0)  '1^  r  (H|,f^) 

2  max  1  h  0  n  max  1 


). 


(A25) 


and 


a  fl_ 

-  ■  ■  ”...  Av  0  ( 


IP9I 


e  r  (H,,n)  '2''  Yj’ftna  r  (il,,n) 

max  i  h  n  n  max  1 
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where 


Gi(y)  =  [  / 


max 


2.  y^2- 
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X  (-?x  In  X  -  v'^) 
min'  '  ' 


iJlMMMMilMM 


IM, 

G2(y)  =  y  Gj(y)  / 


max 


dx 


min 


f  -  2.  2072  * 

x(-2x  In  X  -  y  J 


Here,  ’Snax^y^  solutions  of 

2,  2  2  /.o-,-. 

-  X  In  X  •  y  .  (A27) 

(111)  Density  of  the  group  of  particles  with  the  same  U. 

D 

The  group  of  particles  with  given  Qg  Includes  particles  with  angular 
momentum  Pg  ranging  from  zero  for  orbits  through  the  axis  to  some  maximum 
value  for  circular  orbits.  The  orbit  turning  points  likewise  depend 
on  Pg,  and  In  the  small-orbit  limit  (A15),  we  can  analytically  calculate 
r^x  ®  function  of  flg  and  P0  by  combining  (A20)  and  (A22).  This 
procedure  Indicates  that  all  particles  with  frequency  Og  are  confined  to 
the  region  0  <  r  <  R(fig),  where 

R  ^(£Jg)  -  3a^(l-|ng|/«^)  .  (A28) 

Using  this  result,  (A20)  can  be  rewritten  In  the  form 

I  ?  I  “  I  +  S  +  2(l  +  S  -  2S^)^^^],  (A29) 

TmlnJ 
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0  <  S  =  9|Pg|/(8Y^jMR^QQ)  <  1 

? 


where 


^2 

We  note  that  only  particles  with  S  ■  1/2  and  hence  |Pg|  ■  4y^MR 

pass  through  R,  and  the  orbits  of  these  particles  are  Intermediate 

~2 

between  radial  (Pg  «  0)  and  circular  (|Pgl  “  RY^^MR  All  other 

particles  have  r  <  R. 

max 

In  the  opposite  limit  (A23),  we  can  rewrite  (A26)  as 


r 

max 


(H^.Pg) 


a 

n  0 


'’’e' 


Y.Mll*a  r  (H,,0) 
b  On  max  i’ 


). 


(A3n) 


where 


J  max  I 

with  Xj^jj  again  given  by  (A27).  The  largest  radius  for  a  given  value 
of  flg  Is  found  by  numerically  calculating  the  maximum  f;3(y).  This  gives 

R(ng)  -  1.52  .  (A31) 

For  intermediate  values,  P(^0)  must  be  calculated  numerically, 

and  the  result  for  a  beam  with  a  Bennett  profile  (27)  is  shown  In  Fig. 

2.  We  find  for  all  R  values  that  particle  orbits  with  r  ■=  R(Jia)  are 

indx  0 

Intermediate  between  circular  and  linear  orbits. 

When  Jj^f^(r)  Is  peaked  on-axls,  so  that  Jig  Is  a  monotonlc  function 
of  for  constant  Pg,  the  equilibrium  beam  current  density  (9)  can  be 
rewritten  as 

HgCR-r) 

J^^(r)  -  2  /  dQg  6j^(r;ng),  (A32) 
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where 


6jp(r;a0) 


2qBc  t 
r  i 


max 


dP, 


min 


M(,  I, 

o 


(a:^3) 


Is  the  current  density  of  particles  with  a  ^Iven  precession 
frequency  flg.  In  (A33),  and  P^j^^  are  the  maximum  and 

minimum  Pg  values  of  particles  with  precession  frequencies  equal 
to  (Ig  and  orbits  passing  through  r. 

In  the  small-orbit  limit  (A15),  the  density  profile  can  be 

calculated  In  closed  form.  To  find  the  Pg  limits,  we  use  (A22)  to 
rewrite  the  turning  point  condition  (A12)  In  the  form 


0  . 


(A34) 


We  solve  (A34)  for  |Pgl.  and  using  (A2fi),  we  can  write  the  solutions  as 


1  7  1/7 

"T  V*"  %  ^  - %. - ^  ^  ^  • 


(A3.S) 


Since  |Pg|  must  be  positive,  we  take  »  P^  and  P^^^  ■  max(0,P_), 

noting  that  P^j^  in  (A33)  Is  zero  for 


r^  >  I  R^(ng). 


(A36) 


For  the  Maxwell-Bolt zmann  distribution  (64),  we  also  have 


"  J^p(f'>/(2PY^,MqBcP^)  , 


(A37) 
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2  2 

to  lowest  order  In  r  /a*  The  £J^(r;Q^)  expression  (A33)  can  then  be 

n  do 

written 
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0  max(0,P_)  (P^.- 


(A18) 


which  Rives 


6Jo(r;nQ)  - 


6J, 

6Jo(n;n0) 
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7  ft  ~7 
0<  r^<  ~ 

I  p^(ne>  < 

^2  2 

R  (n.)  <  r^. 


(ASR) 

where  6j^(O;Jl0)  “  ftBJjj^(0)/(3JJ^),  The  profile  (A39)  is  shown  in  PiR. 

1.  The  discontinuity  in  6j^  that  occurs  at  the  outer  radius  R  is  a 
Reneral  feature  that  has  important  consequences  on  the  analytic  structure 
of.  the  elRenmodes. 


Appendix  B!  Evaluation  of  the  Vlasov  Time  Integral 

In  the  Vlasov  perturbed  current  expression  (12),  the  accumulated 
effect  of  the  wave  on  particles  with  some  given  and  is  expressed  in 
terms  of  the  time  Integral 

I(r;H^,Pg)  I  dz'  A(r')  exp[l(e'  "  (*'“  *>]•  ("H 

Here,  z  ■  Bet  is  a  measure  of  propagation  time,  and  r'  and  6'  are  the 
unperturbed  radial  location  and  azimuthal  angle  at  time  z"  <  z  of 
particles  with  the  specified  H^  and  P^  that  reach  the  point  (r,6)  at  time 
z.  The  +  subscript  indicates  the  sign  of  the  radial  velocity  v^  at  time 

z. 

To  put  (Bl)  into  a  more  useable  form,  we  consider  the  change  along 
the  particle  trajectories  of  the  quantities 

A  -1  ^  ^ 

l^(r,z)  -  l[r^(z)]  exp[l^^(z)]  ■  /  dz"  A(r')  exp(l*')  .  (B2) 

Here,  the  phase  angle  is  defined  by 

♦+(z)  -  e_^(x)  - 

and  all  orbit  quantities  are  understood  to  depend  on  Hj^  and  P^,  Since  v^ 
and  6  are  functions  of  r  only  for  these  particles,  we  can  readily 
relate  I_^(r)  to  the  values  at  the  previous  turning  points: 

i^(r)  expd^^)  -  exp(U*^^)  +  f  dl'  expCl^;)  (B4a) 

’^mln 
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(B4b) 


I_(r)  exp(l*_)  »  ^  exp(l^'), 

^max 

where  r^u  and  are  the  orbit  turning  points,  and  we  have  Introduced 

the  compact  notation 


’^mln 


♦♦'".in* 


♦+('*  ■  I 


dr' 


r  ,  "'r*'*’' 
min 


[e(r')  -  fl] 


(B5a) 


*mx  "  '^-^’^max^  ‘  ^  hTTHr  ’ 

r  r 
max 


(B5b) 


and 


dr'  = 


dr' 


|Vj.(r')| 


ecA(r')  . 


(B5c) 


We  eliminate  from  (B4)  by  noting  first  that 


-  Cn  ’  -  ■  ♦-('^mln^i* 


(B6) 


The  relative  phase  between  and  ♦_  Is  fixed  by  taking  the  phase  angles 
to  be  equal  at  r„^„.  so  that  5  We  then  rewrite 

(B6)  as 


♦_(r)  -  ♦^(r)  (B7) 

and  substitute  this  expression  for  In  (B4). 

The  turning  point  values  I  5  ®‘“*  ^mlx  =  ^-^’'max^  **‘® 

A 

determined  by  evaluating  the  full  change  In  between  turning  points: 


I  exp(i^  )  -  I  .  expdi  4_)  ■  /  dr''exp(l^') 
max  max  min  rain  '  + 

min 


(Bfta) 


max 


W*P^‘*mln>  -  ♦mx^^  "  ^  ‘’^'®*P[^^2*mln"  ♦+^1* 


min 


(BRb) 


C  + 


Since  the  continuity  of  particle  orblta  Rlvea  ■  ^nin  ^  ^mln 


*  + 


I  -  1  HI,  (BS)  can  be  solved  for  I  ,  and  I  : 
max  max  max  min  max 


*  T)^exp(-IA^)  +  l>_exp(lA+) 

^raln  "  exp(lA^)  -  exp(-lA^) 


(BSa) 
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P+  + 


max  exp(lA^)  -  exp(-lA^)  * 


(BSb) 


where 


and 


max 


/  dr'exp[+  !(♦'  - 


min 


A*  =  ♦*■  -  4)*'. 

^  max  min 


max 


r"  ivifyT 


»e-“ 


r  j  r 
min 


We  substitute  (B9)  Into  (B4)  to  obtain  expressions  for  l_j_(r)  In 
terms  only  of  orbit  Integrals  from  t^ln  f®  *"  £  ’^max  *  Suppressing  the 
superfluous  +  subscript  on  we  find 
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I^(r)  -  /  dr'exp[l(*'-  ♦)] 
r 

max 

’^max  ,  ^ 

/  dr'{exp[l(4i'-  ♦  -  A4»)]  +  exp[-i(*'+  ♦  -  A^)]}  (BlOa) 
’^min 

A  max 

I_(r)  -  /  dr'exp[-l(^'-  ♦)] 

*^mln 

r 

+  ~2i*^ain( AA)  ^  dT'fexp  [i(<|»'+  ♦)]  +  exp[-l((|>"-  4>)]}.  (BlOb) 

’^mln 

A  1  A  A  ^ 

Since  I(r)  =j  [l^(r)  +  I_(r)]»  we  can  combine  these  expressions  into 

r 

«  ,  max 

Ur)  * /  dr'{exp[i  sgn(r-r')  (♦'-  4>)  ] 

^min 

- [exp(-iA^)cos( t'-  ♦)  +  cos(^‘'+  ^  -  A^)]}  .  (Bll) 

When  the  exponentials  and  cosines  in  (Bll)  are  expanded,  some 
cancellations  occur,  and  we  obtain  the  more  compact  form 
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2i  sin(A^) 


and 


I(r)  ■  -iBc  / 


max 


dr 


^  A 


Iv  (r")l 
r  ,  r 
min 


Trr  A<r') 


X  [  +  0(r-r')8in^co8^‘'  +  0(r‘'-r)8int'co8t]  , 


(B12) 


58 


where  6  is  the  unit  step  function  and  (BSc)  has  been  substituted 

A 

for  dF'’.  This  form  of  I  can  conveniently  be  used  either  to  evaluate 

A 

J(r)  numerically  or  to  Investigate  Its  analytic  properties. 
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